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I. INTRODUCTION 



There has been growing interest recently in studying the relation between gravitational 
field equations describing bulk spacetime dynamics and horizon thermodynamics. A pioneer 
work on this topic was done in jl| where Einstein's equation emerges as an equation of state 
from the basic thermodynamic relation in the Rindler spacetime. A lot of effort has been 
spent to quantify this relation including studies in the so called /(R) gravity |2J and the 
scalar tensor gravity |3] . Moreover, it has been found that the relation between gravity and 
thermodynamics exists also in other spacetimes, including a general static spherically sym- 
metric spacetime [4], the dynamical Vaidya spacetime ^ and cosmological spacetir 
where more modified gravity theories like Lovelock gravity [4] and braneworld gravity |l 0Ml2 1 
have been described by the first law of thermodynamics on the black hole horizon or cosmo- 
logical apparent horizon. Recently, it was further found that for arbitrary diffeomorphism- 
invariant gravity theories, the field equation can be obtained as a state equation of Rindler 
horizon thermodynamics 13l-ll5|. The disclosed relation between horizon thermodynamics 
and bulk gravitational field equation can shed the light on holography [16| and even may 
change the understanding of gravity. In fact, the puzzling thermodynamic feature of gravity 
and/or spacetime is one important motivation of the proposal that gravity might not be a 
fundamental interaction but rather an emergent large scale/numbers phenomenon 

In thermodynamics, the equations governing the equilibrium state of the system can be 
obtained by the extremization of the entropy function. In terms of the connection between 
gravity and thermodynamics, it seems natural that in the gravity context one could also 
define an entropy function(al) for spacetime and derive the gravitational field equation from 



it to describe the spacetime. A novel approach to realize this idea was presented in [18|, |19|, 
where they made an analogy between the spacetime and an elastic solid. The idea of 
spacetime-solid analogy is not new but has a long history, since Sakharov's induced gravity, 
see [20] • At the microscopic level of the analogous picture, its associated degrees of freedom, 
the 'spacetime atoms', are still elusive, but its existence is crucial to interpret the thermal 
phenomenon associated with spacetime. At the macroscopic level, what may be really 
interesting is that, one can expect to develop a theory of elasticity phenomenologically to 
describe the spacetime solid, envisaging gravity as analogous to the elasticity. Obviously, 
the fundamental dynamical variable of spacetime solid is the displacement vector field which 
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describes the elastic deformation of the solid, while not the metric, which may be thought 
of as a coarse grained description of the spacetime at macroscopic scales, somewhat like 
the density of a solid. Although the analogy is not strong, Padmanabhan et al. HQ 
successfully constructed an entropy function of a displacement vector field and obtained 
the correct field equations of Einstein and Lovelock gravity theories by a special variational 
principle on the entropy function. Moreover, when evaluated on-shell for a solution admitting 



a horizon, the extreme value of the entropy agrees with standard Wald entropy [21|, |22] . This 
justifies that the constructed elastic entropy has the significance of gravity entropy. Beside 
further supporting the deep connection between gravity and thermodynamics, this approach 
is attractive since it has the following two ingredients: (a) the metric is not the dynamical 
variable; (b) the field equations remain invariant under the shift T ab — > T ab + Xg ab of the 
matter energy-momentum tensor T ab . It was pointed out that these ingredients may be 
important to solve the cosmological constant problem Jl9|. 



The spacetime-solid analogy formalism proposed in [18|, ll9] leads to the gravitational field 



equations restricted in Einstein theory and Lovelock theory. The standard Wald entropy 
in 18[ is obtained only on the local Rindler frame. Considering that the thermodynamic 
derivation of field equation can be applicable to arbitrary diffeomorphism-invariant gravity 
theories, it is natural to ask: Does the proposed formalism hold in all spacetimes beyond 
Einstein gravity and Lovelock theory? Does it imply something in deep? In this work 
we are going to address these questions. We will analogize the spacetime as an inhomo- 
geneous elastic solid, and develop a general formalism to obtain the theory for arbitrary 
diffeomorphism-invariant gravity. Taking the homogeneous limit, our formalism can reduce 
to obtain the Einstein and Lovelock gravity theories got in 

The organization of the paper is the following. In section II, we define the entropy function 
of a displacement vector field in the spacetime solid for arbitrary diffeomorphism-invariant 
gravity. Then in section III, we extremize the entropy function to derive the equilibrium 
equation. In section IV, we compute the extreme of entropy function on general static and 
non-static spacetimes and show that under appropriate circumstances it is identical to the 
expression for the Wald horizon entropy. We conclude our results in the last section. 
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II. ENTROPY FUNCTION FOR ANY DIFFEOMORPHISM-INVARIANT GRAV- 
ITY 



Now we will describe the spacetime as a deformed solid at the macroscopic level. Usually, 
to describe a general deformable solid with the theory of elasticity, one can define the 
thermodynamical function, like the entropy (or free energy, internal energy etc.) with the 
displacement vector field £ a (x) which describes the elastic displacement of the solid through 
the equation x a — > x a + ^ a (x), to capture the relevant dynamics in the long- wavelength limit. 
Varying the thermodynamical function with respect to £ a , one can obtain the equilibrium 



equation of £ a (see the standard book of elastic mechanics by Landau and Lifshitz 23|). 
Analogizing spacetime to solid, the equilibrium equation of spacetime solid can be expected. 
However, what is remarkable is that the field equations will appear from the equilibrium 



equation 18] . Before achieving this, the key task is how to specify the thermodynamical 
function of spacetime. 

There are several restrictions for the form of thermodynamical function. At first, the ther- 
modynamical function of the spacetime solid should be a scalar, preserving the covariance 
of equilibrium equation. Second, in the theory of elasticity, the thermodynamical function 
can be written as an integral over a quadratic function of small strain tensor, preserving the 



function density to be translationally invariant. In 18] the entropy function of spacetime 
solid was proposed as 

S g [£] ~ f d D x^~g {P abcd V £ b V d i b ) , (1) 
Jv 

by treating V a £& as the strain tensor and P abcd as the elasticity modulus. However, in the 
presence of non-gravitational matter distribution in spacetime, one can not demand the 
translational invariance of entropy density. Hence, the entropy density can have quadratic 
terms in both the derivatives V a £& as well as itself. So it was proposed that the total 
entropy function should also include the matter entropy function 18J] 



S m [£] ~ / d D xV=gT ab U 



b- 

V 



Assuming the "constancy" conditions for elasticity modulus V ' d p abcd = analogous to a 
homogeneous (and isotropy) solid, it was argued that the total entropy function can be used 



to describe the Einstein and Lovelock gravity 



18| 



In our study here we are going to release the constraint V^-P = 0, which is analogous to 



treat the spacetime as an inhomogeneous solid. By this nontrivial extension we will achieve 
to obtain gravitational field equations for more general gravity theories. Before constructing 
the concrete form of entropy function, it is natural to expect to get some hints from the inho- 
mogeneous elastic theory. In mechanics of materials, the single and homogeneous materials 
can be fabricated as piezocomposite materials. The interface between materials produces an 
uneven distribution of stresses which reduces the electric-field-induced displacement char- 



acteristics. 
Materials 



t is interesting to see that these materials, the called Functionally Graded 



24J, have a total free energy: 

F = j^x riptf% ejH - e^ k e l3 E k - ^E t E^j , (2) 

where denotes the components of the strain tensor, Ei is the component of the electric field 
vector related with the displacement vector by the constitutive equations. The quantities 
pijki anc j represent the elastic and the piezoelectric modulus, respectively, and is 
the electric permittivity. Motivated by this free energy, we present that a general entropy 
function of spacetime solid should be 

S g [£] ~ / d D xy/=g (P abcd V c ^d^ b + V d P o6cd V c £ a 6 + V c V d P acbd U b ) . (3) 
Jv 

Thus roughly one may make the analogy of the last two terms in (J3]) as the piezoelectric 
effect of inhomogeneous spacetime solid. 

We need to determine the elasticity modulus p abcd ^ especially its symmetry which, in the 
usual elastic solid, will characterize symmetry of structures of crystals, such as monoclinic 
and tetragonal, etc. However, the concrete form of P abcd only can be determined in a 
complete theory by the long wavelength limit of the microscopic theory just as the elastic 
constants can in principle be determined from the microscopic theory of the lattice. In 
macroscopic level, one can only know that the structure of the gravitational sector is encoded 
in the form of p abcd y so the object P abcd should be built out of metric and other geometric 
quantities. Such a tensor can be constructed as a series in the powers of the derivatives of the 
metric. One may expect that the lowest order term leads to Einstein's theory while hi ghe r 
order terms come from the quantum corrections of underlying microscopic theory. In 18], 
assuming p abcd has symmetry with Riemann tensor R abcd and impose V ' d p abcd = 0, the lowest 
and second order correspond to Einstein's theory and Gauss-Bonnet theory (then Lovelock 
theory), respectively. However, if we release the condition V ' d p abcd = 0, more general 



tensor are possible. The general construction of P abcd should be universal for arbitrary 
diffeomorphism-invariant gravity, including the tensor for Einstein gravity and Lovelock 
gravity as our special cases. Assuming the symmetry of P abcd in identical with R abcd , a 
simplest candidate is 

P abcd = 5=-, (4) 

Oriabcd 

where L is the Lagrangian of gravity theories. 
Finally, the total entropy can be written as 

S[i] = S g [£] + S m fc] (5) 
= 4 J d D xV=g (V bcd v c £ a V d & + V d P abcd V c Ub + V c V d P acbd U b - ~T a6 £ a £^ , 

where some proportional constants are chosen with hindsight. It should be pointed out 
that the entropy function Eq. (TjQ) for Einstein and Lovelock gravity and their generaliza- 
tion Eq. (jSJ) for general gravity can not be fixed completely by the analogy with elastic 
theory. One can understand that the entropy functions Eqs. flTJ and flSJ are constructed 
phenomenologically. One of the desired phenomenon is to obtain the field equation from the 
equilibrium equation. The other is to justify the entropy function with the significance of 
gravity entropy, which is implemented by identifying the on-shell entropy function with the 
Wald horizon entropy. We will show both of them in sections below. 



III. FIELD EQUATIONS FROM EXTREMIZING THE ENTROPY 

In this section, we will derive the equilibrium equation of spacetime solid. Although the 
expression (jSj) is well defined for any displacement vector field, one can only obtain significant 
results for suitably chosen vector fields. The vector is required to characterize the special 
property of spacetime solid. The most nontrivial property of spacetime is the existence of 
the horizons which act as one-way membranes which block information for a specific class 
of observers. The existence of horizons, which are null hypersurfaces, is a feature of any 
geometrical theory of gravity and is reasonably independent of the field equations. We 
hence assume that the spacetime solid is deformed induced by the change of horizon. As 
the simplest case, one can consider that the matter is freely falling into the horizon along 
the transverse invariant ingoing geodesies. Alternatively, one can also consider a virtual 
displacement of horizon radially normal to itself engulfing the matter. Obviously in this 
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case the total displacement of spacetime solid is induced by the change of horizon. On the 
true horizon, this change, i.e. the displacement vector £ , should be characterized by its 
outward null normal vectors. Actually, we will derive the correct field equations when we 
impose £ a as null vectors and obtain the entropy of horizon when we use the outward unit 
normal of a near horizon surface to approach the null normal of true horizon. 

Varying the vector field £ a after adding a Lagrangian multiplier X(x) for imposing £ a with 
constant null norm 5 (£ a £ a ) = 0, we find: 

6S[£] = 4 / d D x^[2P abcd V c ^VMb + V d P abcd V c £ a <5& + V d P abcd V c 5Ub 
Jv 

+V c V d P acbd (5Ub + - \ {T ab + Xg ab ) 

= 4 [ d D x^[2V d (P abcd V c ^ b ) - 2P abcd V d V c ^ b - 2V d P abcd V c ^ b + V d P abcd V c ^ b 
Jv 

+V c (V d P abcd 5Ub) - V c V d P abcd 5Ub - V d P abcd 5i a V c i b - V C V, (2P acdb + P abcd ) 

- l -{T ab + \g ab )iMb]. 

where we have used the symmetry of Riemann tensor p abcd = pNN, pabcd _ pcdab^ ^ n( j 
we get the second to the last term above following 

V c V d P acbd (5Ub + = - V c V d (P cabd + P cbad ) i a 5i b = - V C V, (2P acdb + P abcd ) £ a 8£ h , 

(6) 

by using rotational symmetry of p a \- bcd i = 0. Then 6S[£] can be simplified as 

6S[£] = 4 / d D xV=g[2V d (P abcd V c Z a 5£ b ) + V c (V d P abcd 6Ub) - 2P abcd V ' d V ' c ^ b 
Jv 

-2V c V d P acdb U5i b - l - (T ab + \g ab ) 
which leads, under the stokes theorem, 

6S[£] = 4 / d^xVhpn^P^VcZaS^) + n c (V d P abcd 5Ub)} 

JdV 

-8 J d D x^[P abcd V d V c iJi b + V c V d P acdb i a 5i b + J (T ab + \g ab ) 

where n a is the vector outward normal to boundary dV and h is the determinant of the 
intrinsic metric on dV . As usual, we set the variation S^ a to zero at boundary (Even though 
the £ a is not fixed on the boundary, we still do not care about the boundary term since it 
only gives the boundary condition of £ , which will not affect the bulk equilibrium equation 
which is independent of £ a as we will show.). Therefore the first integration in the above 
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equation vanishes, and the condition that S[£\ be an extremum for arbitrary variations of 
£ a leads 



6S[£] = -8 f d D Xy/=g~ 
Jv 

= 4 / d D x^ r g 



pabcd VdVc ^ b + S7 c V d P acd %6£ b + i (T ab + \g ab ) £ a 5£ b 



P bedc R a edc - 2V c V d P acdb - - (T ab + Xg ab ) 



Jv 

where we have used the definition of the Riemann tensor in terms of commutator of covariant 



Za5£ b , (7) 



derivatives, as well as some alteration of index. One can find that the equilibrium equation 
that follows from our variation principle of entropy function is 

P bedc R a edc - 2V c V d P acdb - ~T a6 - Ixg**) £ a = 0. (8) 

Requiring the condition (jSJ) to hold for arbitrary null vector field (It is interesting to note 
;hat this requirement is also invoked in the derivation of field equation as a state equation in 
l|), one finds that the equilibrium equation from extremizing the total entropy is reduced 



to 



pbedc R a edc ~ 2V c V d P acdb - ^T ab - ^Xg ab = 0. (9) 



This is a remarkable result that we have obtained a dynamical equation governing the 
background instead of the usual equilibrium equation determine the displacement vector 
field. This uninstinctive situation happens because the symmetry of tensor p abcd and the 
entropy function are so special that Eq. (JSJ does not contain derivatives with respect to £ a . 

By demanding conservation of the stress tensor and using the Bianchi identities, one can 
find that the Lagrangian multiplier 

X - L ~^G' < 10 > 
where A is an integration constant and G is the Newton gravitational constant. One can 
immediately see that the equation ([9]) is just the exact field equation derived from ordinary 
variational principle from the action of arbitrary diffeomorphism invariant theories of gravity 



22| (involving no more than the second derivatives of the spacetime metric g a b and the first 



derivatives of the matter fields \&„ 



/ d D xy^gL(g ab ,R abcd ,^ m ,V a ^ m ), 
Jv 



supplemented by appropriate generalizations of Gibbons-Hawking-like boundary terms. One 
can find two crucial features of our derivation that, the variational principle is based on 
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a vector with constant norm in spacetime instead of the usual metric field g a b, and the 
equilibrium equation is invariant under the shift T a b — > T a b + Xg a b since Xg a bC, a ^ b = Xe is not 
varied when £ a is varied, regardless of whether \7 d p abcd = Q. 



IV. ON-SHELL ENTROPY FUNCTION 



The result in the previous section provides an alternative variational principle in deriving 
field equations of arbitrary diffeomorphism invariant theories of gravity. In the following, we 
will show that the boundary term of entropy function S \£] will lead to the standard Wald 
entropy. The specific case on this topic was discussed in jl8l ] . 

Manipulating the covariant derivatives of Eq. (]5]), the entropy function can be rewritten 

as 

S[€\ = ^j v d D x^g(P abcd v c £ a v d tb + v d P ahcd v c u b + v c v d P acbd Ub - \r ab u b ) 

= 4 ^ d D x^V d (P abcd \7 c Ub) ~ 4 J d D x^g~(P abcd V d V c Ub + V c V d P acbd Ub - -T ab Ub) 
= A f d D ~ l xn d (P abcd V c Ub) + 2 / d D x^g-(P bedc R a edc - 2V c V d P acdb - ]-T ab )^a 

JdV JV * 

= 4 / d D ~ 1 xn d (P abcd V c Ub)+ [ d D x^gXe, (11) 

JdV JV 

where similar derivations of Eq. OH]) has been used in the second equality, and the field equa- 
tion (Q has been used in getting the last line. This result shows that, as mentioned in 15), 
the second term of third line can be thought of as the entropy in any diffeomorphism- invariant 
gravity theory, except a total divergence. Then one can recover the entropy function Eq. 
((5]) by reversing the derivation in Eq. ( TTTT) . 

One can find that the final result of Eq. (ITT]) is not affected by \/ d p abcd ^ o. The second 
term of the last line vanishes since e = (Even if e is a nonvanishing constant, it is not 



a surface term 



18 



19)). Therefore, we will concentrate on the first term, which will be 
interpreted as the surface entropy of horizon. 

At this stage, we have not put any restriction of the boundary dV . In general, the 
boundary is (D — l)-dimensional. Hence the entropy function should not be the desired 
entropy of some (D — 2)-dimensional section of horizons. However, it was pointed out in 
[18| that, when part of the boundary dV is null, it is intrinsically (D — 2)-dimensional. This 
case needs to be handled by a limiting procedure from near horizon to true horizon. 



In [18J , the standard Wald entropy of Einstein and Lovelock gravities are obtained from 
the boundary term on the local Rindler frame through a limiting process to the null Rindler 
horizon. They also mentioned that the same result can be recovered for any static spherically 
symmetric spacetime. Here we will extend their discussion to the general static spacetime 
without spherical symmetry and even to the stationary but non-static spacetime. 

We will briefly introduce the coordinate system which is suited for the discussion of the 



general static spacetime (see the details in 25j). Given the general static spacetime, one 



can decompose the metric into a block-diagonal form as follows 

ds 2 = —N 2 dt 2 + g^dx^dx" , fi, v = 1, 2, • • • 

We can arbitrarily choose a particular (D — 2)-surface in the constant-time slice and utilize 
Gaussian normal coordinates in the surrounding region, 

g^dx^dx 1 ' = dn 2 + gABdy A dy B , A, B — 2, 3, • • • 

where n represents the spatial direction normal to the specified (D — 2)-surface. The Killing 
horizon H, generated by the timelike Killing vector field x = dt is approached as N 2 — > 0. 
One can verify that k = d n N\ Q complies with the standard version of the surface gravity. 
This enables us to write a near horizon Taylor expansion for the lapse N(n, y) = nn + 0(n 2 ). 
Since we want the horizon to be regular and not possess a curvature singularity, some 
curvature invariants must remain finite in the horizon limit, which enables us to refine the 
expansion for the lapse as 

N(n, y) = nn+ ^-n 3 + 0(n 4 ), (12) 

and write 

9A B (n, y) = [g H ] AB (y) + M^M^ 2 + 0(n 3 ). (13) 

One can find that the metric (Tl2|) and (fl3|) are more general than the metric of local Rindler 
frame in the presence of K2 and [g<2\ ab . The later directly leads to the nonvanishing component 



of Vb£d, not only V*£t> which makes the identification in [18| become invalid between the 
entropy function and Wald entropy through introducing binormal to the cross section of H. 

We shall evaluate the surface integral for S [C] on.— sfeeii near Killing horizon if on a timelike 
surface S, which is denoted as n =constant and the neighboring spacetime is describe d by 



metric ffl2|) and ffTB"]) . This timelike surface £ can be called as the stretched horizon 
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26 



23, 



which has a non-singular induced metric and then provides a more tractable boundary than 
the true horizon. A rigorous one-to-one correspondence between points on the true and 
stretched horizons can be realized by using, for example, the null rays that pierce both 
surfaces. Roughly, we can expect that the on-shell entropy will match the Wald entropy 
under the limit n — > in the end of the calculation. Take £ a = n a as the spacelike normal 
to these surfaces E, with components 

^ = n a = (0,1,0,0,...) (14) 

and unit norm. In the limiting process the spacelike unit vector will be a null vector, which 
denotes that we are considering the null surface as a limit of a sequence of timelike surfaces. 
The metric determinant h of these surfaces E can be decomposed as y/h = Nyfa, where a 
is the metric determinant on the transverse spatial surfaces, having the limit on the true 
horizon yfcj — > yfgu. 

Now we will check whether the surface term 

S Pad = 4 / d D - l xVhn d (P abcd V c Ub) (15) 
will be reduced to Wald entropy in the limiting process. At first, we will introduce the Wald 



entropy |2l|, |22] based on a simplified version of the formalism [28]. Consider a generally 
covariant Lagrangian L, that involves no more than quadratic derivatives of the spacetime 
metric g a b- Under the diffeomorphism x a — > x a + \ a the metric changes via Sg a b = —V a Xb — 
V a Xb- By diffeomorphism-invariance, the change in the action, when evaluated on-shell, is 
given only by a surface term. This leads to a conservation law, V a J a = 0, for which we 
can write J a = VbJ ab - Here J ab defines (not uniquely) the antisymmetric Noether potential 
associated with the diffeomorphism x a , which can be formulated as 

J ab = -32tt {P abcd V cXd - 2 Xd V c P abcd ) . 

Associated with a rigid diffeomorphism x a , there is the Noether charge defined by integrating 
the Noether potential over any closed spacelike surface S of codimension two. It turns out 
that the corresponding Noether charge is just proportional to the entropy 

1 ' Tab, 



Swaid — t: / J a dB ab , (16) 
8k Jb 

when x a is a timelike Killing vector (with vanishing norm at the Killing horizon), and the 
surface B is the cross-section of Killing horizon H. However one can formally define the 
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quantity Swaid on any closed spacelike surface f2, and only in the end take the limit in which 
Q approaches a section B of the Killing horizon H. In the following, we will define such a 
quantity on the section Q of a stretched horizon S described by metric (1121) and (ITS]) , and 
compare Sp a d with Swaid i n the limit n — > 0. On the stretched horizon £, the timelike 
Killing vector can be specified as 

X a = (1,0,0,0,...)- (17) 

The proper velocity u a of a fiducial observer moving along the orbit of x a is u<1 = 
X° ' I \/—gabX a X b = (i^) a wnere r is the proper time. The fiducial proper velocity u a and 
unit normal n a of stretched horizon £ define dVL ab = n^ a u )dA, where dA = ^d D ~ 2 y is the 
area element on cross section Q. We will evaluate both Sp arf and Swaid for several typical 
diffeomorphism invariant theories of gravity. For simplicity, we will restrict on D = 4, but 
the results can be directly generalized to more higher dimensions. 



A. L~f(<j>,R) 



As the first example let us consider the following Lagrangian 

1 



16nG 



Obviously, this example contains the popular f(R) and scalar-tensor gravity as its special 
cases. From the Lagrangian the tensor p abcd reads as 



->abcd 



dL 



1 df 



{g ac g bd - g ad g bc ) 



(18) 



dR abcd 327rGdR 

Substituting this tensor and the normal vector Eq. (|T4l) into Eq. (Tl5l) . and preserving the 
leading term of n in the end of the calculation, we obtain 

Spad = 4 / d 3 xVhn d (P abcd V c n a n b ) 

+ O {n 2 ) 



,3 ,— 9f 



8ttG 



Restricting the t integral within the range (0, 2ir/n) for the periodicity in Euclidean time 



33 



341 ]. we can obtain 



S 



Pad 



_l_dl 

AG OR 



gnd 2 y. 
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When /(</>, R) = R, the entropy density is -^q, and Sp a d is the well-known Gibbons- 

Hawking entropy Substituting Eq. (|T8l) . Killing vector (|T71) . and dQ ab = n( a u )dA into Wald 
entropy (TL6I) . and preserving the leading term of n at last, we get 

(19) 



Swald — T— I J^dQab 



K Jn 
udR 



We notice that although the higher-order terms O (n 2 ) of Sp a d and Sy/aid are different, their 
leading terms are exactly the same. 



B. L ~ aR abcd R abcd + PR ab R 



jab 



Beside the term f(R) ~ R 2 , the higher (quadratic) derivative interactions usually include 

(o:RabcdR abCd + /3RabR ab ) , 



1671G 

with arbitrary parameters a, (3, which may depend on some scalar fields. Its derivative with 
respect to R a bcd is 



jabcd 



8nG 



aR abcd + _p ( g bd R ac _ g ad R bc + g ac R bd _ g bc R ad^ 



Similar to the above case, we can obtain 



S 



Pad 



, )r , ([92] 22 [g H ] 33 + [g H ] 2 2 [#2] 33 + [g H } 22 [^] 33 - 2 [g H } 23 [^] 23 - MS?) 

n ^gH k k 

2a^2 



+ 0(n 2 )}^d 2 y 



Gk 

which is identical with the evaluation of Sy/aid up to O (n 2 ). 



In 



18], it was proved that the Wald entropy of Gauss-Bonnet gravity L ~ R abc( iR 



abed 



4R ab R ab + R 2 is identical with the entropy derived from Eq. ( fTBT) (in Rindler frame). We 
have checked that this is a special case of combination of our discussions in A and B, since 
the formula Swaid and Sp a d are both linear in the Lagrangian. 



C. L ~ R ab V a ctN h ct> 

It is known that there are some ambiguities in Wald entropy. Considering the followin 
interaction involving the metric and a scalar field 
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g 



the corresponding p abcd reads 

pabcd = ^v [ W ][c VV- 

167TG 

It is important to notice that this tensor is not altered if one adds some terms about scalar 
fields (but no more than their two order derivative) into L, such as, 

1 16nG 

In 



v a vVv a v„0 - (v 2 <py + /uv>v 6 

291 ]. it was pointed out that, L^e (e is the volume form) can be written as a total derivative 
L^e = dcti, which yields a vanishing contribution to the Wald entropy. This contradicts with 
the direct evaluation using the tensor p abcd in Wald entropy on the section of Killing horizon, 
until one realizes that L x = on it. For the entropy Sp a d, one can evaluate 

S Pad = --^ / d 3 xVhn d (V^g b ^ c V d ^V c n a n b ) 

| + 0(n) V9^d\ (21) 

which is identical with the evaluation of Swaid up to O (n), and seems nonvanishing. In fact, 
one can make the similar understanding to that in [29], where S Pad ~ n n W a (t) should vanish, 
since n a will be a null Killing vector when n — > 0, and hence n a V a ~ L x <p vanishes on the 
Killing horizon. 

Hereto, we have shown that the surface term Sp a d can be reduced to Wald entropy near 
the static horizon in the leading order (the higher order terms are different). We will further 
show that this result can be generalized to any stationary but not static black holes. For 



such a black hole, it is expected to be axially symmetric |30j. Consider our spacetime is 
invariant under "time-reversal" , it is convenient to write the spacetime metric near horizon 
n = as 

ds 2 = —N (n, z) 2 dt 2 + (n, z) [dip — u (n, z) dt] 2 + dn 2 + g zz dz 2 , 

where iV denotes the usual lapse function and u is the angular-rotation parameter. The ze- 
roth law of black hole mechanics and rigidity theorem [32[ for axially symmetric (stationary, 
non-static) Killing horizons tell that surface gravity n = d n N\ n ^ and ujh = w| n _>o mus t be 
a non-negative constant on the horizon, respectively. Moreover, a stationary Killing horizon 
is a geodesic submanifold which implies that the horizon is extrinsically flat, i.e. the extrin- 
sic curvature of the horizon must be zero KyJ\ = 0. The above properties directly imply 
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a set of necessary constraints 



g w (n, z) = [g H ] w (z) + !^n 2 + 0(n 3 ) 
g zz (n, z) = [gn] zz (z) + [ ^M n 2 + 0{n% 
u(n, z) = u H + ^±n 2 + 0(n 3 ), 
N(n, z) = nn + Mfl n 3 + 0(n 4 ), 



(22) 



where the first-order terms in uj and TV are required to vanish to avoid a curvature singularity 
on the horizon. The Killing vector of this spacetime is 



X" 



1,wh,0,0), 



(23) 



in terms of (t, f, n, z) coordinate system. 

Invoking the time integration carried out in Euclidean sector, and using the metric fT22l) , 
the Killing vector (|23|) and the corresponding fiducial velocity u a , we can evaluate the surface 
term Sp a d and Wald entropy near the horizon for different gravity theories. We find the 
entropy on horizon: 



Spad — I 



df_ 

ndR 



gndzdip = Swaid + O (n 2 ) , for case A, 



S 



Pad 



^at^T^vp ^» - i \Sh\\ v [9h} zz u 2 {z) 2 + [g H ] [g 2 ] zz + ^ [g H ] w [g H ] 



a 



( 3 [fl'tflL w 2 - 4kk 2 ) + O (n 2 )}y/g^dzdip 



4 AGk 2 
Swaid + O (n 2 ) , for case B 



S 



Pad 



g H dzdip = Sy/aid + O (n) , for case C. 



Now we can conclude that the boundary term Sp ac [ are the same as the standard Wald 
entropy for general static and stationary but non-static black holes. This is one of the key 
results of our paper, which justifies our function S [£] as an authentic 'entropy'. 



V. CONCLUSION AND DISCUSSION 



In thispaper, we have generalized the analogy between the spacetime and solid devel- 
oped in 18). We have shown that the spacetime with generalized gravity theory can be 
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analogized to the inhomogeneous elastic solid. Spacetimes described by the Einstein and 
Lovelock gravity are the limiting cases when we consider the solid being homogeneous. One 
key point in this work is that we have successfully constructed an entropy function of a 
displacement vector field to describe phenomenologically the macroscopic level of spacetime 
solid. Extremizing the total entropy function with respect to the displacement vector, we 
have found that the equilibrium equation can be identified with the field equation of gen- 
eral diffeomorphism-invariant gravity. Our generalization supports that there is an analogy 
between spacetime and solid and that the gravitational thermodynamics is not restricted on 
the concrete gravity theory [13(. We also expect that this approach is helpful to obtain the 
: ield equation in the scenarios of emergent gravity where the correct field equation is absent 



Beside providing an alternative approach to obtain the arbitrary gravitational field equa- 
tion, we have shown that the entropy function on the boundary of any stationary spacetime 
is identical with the standard Wald horizon entropy in several typical higher derivative the- 
ories of gravity, though a general relationship between these two entropy expressions still 
needs further understanding. This provides a new method to calculate the horizon entropy 
for any stationary spacetime besides the Hamiltonian method, Noether charge method and 



the field redefinition method 



29|. 



It should be stressed that the spacetime solid is very special. This solid acquires the 
equilibrium equation independent of the displacement vector field. We suspect that this 
acquirement may be related to certain underlying first principle to constrain the macroscopic 
description of spacetime solid. Moreover, the displacement vector field has the zero norm, 
which characterizes the key role taken by the null hypersurfaces in the spacetime solid 



35 



36|. Furthermore, the spacetime solid has a special surface, the stretched horizon. 



Technically, we have only treated it as the near horizon to approach the true horizon. 
However, in the membrane paradigm, it has been treated as a dynamic fluid membrane, 
obeying such pre-relativistic equations as Ohm's law and the Navier-Stokes equation Q- 
Obviously, it is valuable to study the thermodynamical property of stretched horizon in the 
spacetime solid, comparing the differences and similarities with the fluid membrane. If one 
can understand these nontrivial properties of spacetime solid, the analogy might be of use 
in the context of the semiclassical limit and even the quantum gravity (spacetime) enigma. 
In fact, recently it have been proposed that if the entropy function for Einstein and 
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Lovelock gravity is interpreted as an action in the semiclassical limit, its value will affect the 
phase of the semiclassical wave function, then the observer independent of the semiclassical 
gravity requires this phase (i.e. the Wald horizon entropy) to be quantized in units of 2n. 
Our work further suggests that the Wald entropy can be similarly quantized in arbitrary 
diffeomorphism-invariant theories of gravity. 
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